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Abstract 

Ground state energies and on-site density-density correlations are calculated 
for the 1-D Hubbard model using a linear combination of the Hubbard pro- 
jection operators. The mean-field coefficients in the resulting linearized Equa- 
tions of Motion (EOM) depend on both one-particle static expectation values 
as well as static two-particle correlations. To test the model, the one particle 
expectation values are determined self-consistently while using Lanczos de- 
termined values for the two particle correlation terms. Ground state energies 
and on-site density-density correlations are then compared as a function of U 
to the corresponding Lanczos values on a 12 site Hubbard chain for 1/2 and 
5/12 fillings. To further demonstrate the validity of the technique, the static 
correlation functions are also calculated using a similar EOM approach, which 
ignores the effective vertex corrections for this problem, and compares those 
results as well for a 1/2 filled chain. These results show marked improvement 
over standard mean-field techniques. 
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Following the discovery of high T c superconductivity in the copper oxides [1], strongly 
correlated systems have generated a great deal of interest in condensed-matter theory. It is 
recognized that the electronic properties of many systems, e.g. the heavy- Fermion systems 
and high T c superconductors, are basically determined by strong correlations [2]. The funda- 
mental model describing a tight binding system with strong local correlations is the Hubbard 
model [3], characterized by the hopping matrix element, t, and the on-site interaction energy, 
U. Despite the simple nature of the model (basically a one parameter model, U/t), it is in 
general not analytically tractable. Interestingly enough the model has been solved exactly 
for the 1/2-filled band in one dimension [4]. Solutions have also been obtained under certain 
limiting conditions. In the limit of large U, allowing for only single occupancy at each site, 
expansions in 1/iV, where N is the spin degeneracy have been solved [5] as well as an exact 
mean-field solution in the limit of large dimensions [6]. However, in spite of the physical 
insight gained from these approaches, these 'exact' solutions are clearly approximations to 
the actual physical model which involves a spin degeneracy of 2 and dimensionality ranging 
from 1 to 3. Thus, a full understanding of the physics described by the Hubbard model has 
so far escaped theoretical physicists. What we wish to describe in this communication is 
a straightforward mean-field scheme that linearizes the Equations of Motion (EOM) that 
result from the Hubbard model and then compare the resulting fits to exact numerical Lanc- 
zos solutions for a 12 site system. We start with the charge conjugation symmetric version 
of the Hubbard model, 

H = - tijclcjv + ^ - \)( n i-° - \)- (!) 

The EOM obtained from this model using finite temperature formalism are well known and 
are given by 
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The standard mean-field approximations (which were also employed by Hubbard [3]) are 
n-i- a Cj a ~ (ni- a )cj a , and (c]_ a Cj- a — c^_ a Ci- a *)Ci a ~ (c]_ a Cj- a — c^_ a Ci- a )ci a . Due to trans- 
lational invariance of the system, the second or off-site term is ignored and the mean field 
value of (ni_<j) is given by n e /2 where n e represents the electronic filling of the system. 
The resulting mean-field solution demonstrates the short range Coulomb-generated energy 
gap and the corresponding splitting into two bands but is well known to exhibit several 
shortcomings, the most obvious being the absence of a Mott transition for all dimensions. 

We have decided to expand the usual mean-field linearization by projecting out ni_ a Ci a 
terms as well as terms proportional to c 8(7 . To do this, we form a linearly independent 
combination of the Hubbard projection operators first described by one of the authors [7]. 
Our projection operators then take on the form, 
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This definition of the projection operator satisfies Pc 8(7 = c 8(7 and Prii- a Ci a = rii_ a Ci ai and 
thus is consistent with the idempotent property of projection operators, P 2 = P. Addition- 
ally, it incorporates terms in c 8(7 and rii- a Ci a that the original mean-field scheme ignored. 
Equation 2a above is unchanged but the new mean-field approximation for 2b now reads, 
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where, for a translationally invariant system, 
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Using the standard definitions for the Green's functions in this problem, gij a {j) = 
-(T t c 8(7 (t)c] (7 (0)), and / 4 >(t) = - (T T n t _ a (T)c ta (T)c] a (0)) , the mean-field EOM for the 
Green's functions are conveniently expressed as 
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where Hf^ F is given by 
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The coefficient for the on-site term, 8, can be expressed as simple linear combination of a 
single particle expectation values for / and its equivalent charge-conjugate expression for 
holes. However, the difficulty in obtaining mean-field solutions to the Hubbard model us- 
ing this approach is in obtaining values for the static two particle correlations that occur 
as coefficients for the off-site terms in the linearization process. These correlation coef- 
ficients demonstrate the direct coupling that occurs between the mean-field solution and 
the multiparticle correlations that arise in strongly correlated systems. These correlation 
functions can be calculated using a standard linear response formalism by taking variations 
of the EOM with respect to a density field, yielding density-density correlations, or to a 
bond field, yielding bond-bond correlations. The resulting two-particle correlations are then 
coupled to three-particle correlations. Rather than extend this heirarchy, we have decided 
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to demonstrate the significant improvement of this technique, as compared to the original 
mean-field approach, by using numerically obtained Lanczos wave functions to calculate the 
two-particle correlations. With the coefficients for the off-site terms fixed, this mean-field 
problem can then be solved self-consistently. The solutions for both g and / are then used 
to predict ground state energy, E g , and the on-site density-density correlation, (n,vni_ ff ), for 
a 12 site chain. These values, as well as those obtained for the original Hubbard solution, 
are shown in Figs. 1(a) and 1(b) as a function of U jt for a 1/2 filled system and in Figs. 2(a) 
and 2(b) for a 5/12 filled system. The significant improvement of this mean-field scheme is 
especially noticeable for the 1/2 filled case. The difference in the ground state energy be- 
tween the "Hubbard mean-field" (HMF) and and the Lanczos (L) solution is approximately 
1/2 the difference between the standard mean-field (MF) and the Lanczos result (albeit 
with a different sign). The improvement in the on-site correlations is even more significant 
with the difference between the HMF solution and the exact result being approximately 
1/3 the difference between the MF solution and the Lanczos result. Since the MF solutions 
themselves for 5/12 filling are significantly improved as compared to the 1/2 filled case, the 
observed improvements of the HMF solutions for this filling are not as impressive but are 
none-the-less significant, particularly for the on-site density correlations. 

To further validate the improvement in this mean-field scheme, we can again use the 
EOM approach to calculate the static correlation functions that were used as coefficients in 
this mean-field approach and compare them as well to the Lanczos results. This approach 
ignores all vertex corrections and consequently the three particle correlations that we alluded 
to above. However, for the sake of simplicity, it can further demonstrate the effectiveness of 
this approach. First, we define two pairs of coupled time-dependent correlation functions; 
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as well as 
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Taking r derivatives with respect to r 2 of the above correlation functions leads to sets of 
EOM very similar to those for the Green's functions, 
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Inverting these EOM to obtain the time-dependent response or correlation functions 
yields a set of expressions from which we can evaluate the static correlation functions required 
for our Hubbard mean-field scheme. However, before we proceed it is interesting to note 
some of the self-consistencies that occur with these EOM, most notable among them include; 
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Thus diagonal components of the correlation functions reduce to the correct Green's func- 
tions in this EOM scheme. Meanwhile, the density-density correlation that we require can 
be written, 

(6n l+l _ !T 6n l _ !T ) = Q l+ll+lll -a(00~) - (n_ !J ) 2 = 
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while the bond-bond correlation of interest reduces to, 

( C i+la C iv ( C i-a C i+l-v C i+l-a C i-v )) = 
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G l+lll+lia (00 ) + G l+llll+la (00 ) = -2(8n i+1 - a 8ni- a ). (10) 

Here ^H^ 1 is the 12 component of the inverse of the mean-field Hamiltonian operator. 
These expressions for the static correlations of interest can easily be evaluated from the 
mean-field solutions. These results can also be compared to the correlations obtained from 
Lanczos wave functions. These results for a half-filled 12 site band are shown in Fig. 3. 
They begin to show some differences, albeit with similar trends, for U on the order of the 
band width. This is probably largely due to the omission of relevant vertex corrections 
which become important at larger [/'s, rather than being due to significant errors in the 
mean-field scheme as presented. In fact, self-consistent mean-field results can be obtained 
using this scheme with the above expressions being used for the correlations. The results 
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are virtually identical to those presented here for U < 4t. The point we wish to make here 
is that this mean-field scheme for the Hubbard model represents a significant improvement 
over standard mean-field approximations and can be easily, as well as accurately, applied 
for U less than or equal to the bandwidth ill < 4t). The authors also feel that by including 
vertex corrections while ignoring 3-particle fluctuations, this mean-field scheme can answer 
many questions for the nearly half-filled Hubbard model, including the Mott transition in a 
3-dimensions. Additionally, this approach can easily be extended to include nearest-neighbor 
repulsion effects arising from Vij. 

The authors acknowledge Andrei Ruckenstein for suggesting this problem and the Office 
of Naval Research for supporting this work. 
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Figure la. On-site charge correlations as a function of U for a 12 site 1/2 filled band. t=l is used in all figures. 
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Figure lb. Ground-state energy as a function of U for a 12 site 1/2 filled band. 
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Figure 3. Density-Density correlations and Bond-Bond (B(s) 
lations as a function of U for a 1/2 filled band. 
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